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Abstract 

Defects introduced in ferromagnetic nanodisks may deeply affect the structure and dynamics of stable 
vortex-like magnetization. Here, analytical techniques are used for studying, among other dynamical aspects, 
how a small cylindrical cavity modify the oscillatory modes of the vortex. For instance, we have realized 
that if the vortex is nucleated out from the hole its gyrotropic frequencies are shifted below. Modifications 
become even more pronounced when the vortex core is partially or completely captured by the hole. In these 
cases, the gyrovector can be partially or completely suppressed, so that the associated frequencies increase 
considerably, say, from some times to several powers. Possible relevance of our results for understanding 
other aspects of vortex dynamics in the presence of cavities and/or structural defects are also discussed. 
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1 Introduction and Motivation 

Ferromagnetic Permalloy nanodisks with lateral dimension (thickness, L) around some dozens of nanometers 
have been fabricated and investigated for their potential applications in a number of magnctoclcctronic mecha- 
nisms. In particular, it has been observed that above the so-called single-domain limit, magnetic vortex states 
appear in these samples PQ [2J [3j, exhibiting a planar- like arrangement of spins outside the core, where a perpen- 
dicular magnetization is observed [H 03 [B]. As long as one could manipulate these states other possibilities would 
emerge. In fact, one way towards this control is obtained by removing some small portions of the magnetic nan- 
odisk, in such a way that the defects (cavities) so created work by attracting and eventually pinning the vortex 
around themselves [7J [HI GO EH EQ E^l EH E4] (similar effect also applies to soliton-like spin excitations [T5]). 
Based upon such an idea, Rahm and coworkers [TB] have studied the cases of two, three and four holes (each of 
them with diameter ~ 85 nm) inserted in a disk with diameter ~ 500 nm, separated by around 150 nm — 200 nm. 
Their experimental results confirmed the previous statement about vortex pinning and put forward the possi- 
bility of using these stable states as serious candidates for magnetic memory and logical applications as long as 
we could control vortex position, for example, applying a suitable external magnetic field which should shift the 
vortex center from one defect to another, and vice-versa. Basic logical operations have been obtained in these 
lines by means of bistable magnetic switching[17J. 



When the vortex experiences a suitable alternating field its gyrotropic mode may be resonantly excited in the 
subGHz range, which has attracted a great deal of efforts for it is the lowest translational vortex mode [T8"l[H)ll2T)I 
|2"T1 |2"21 |2"51 l2~i] . Such a frequency depends upon the geometry of the nanodot and, in the case of thin cylindrical 
samples, it was initially believed that the aspect ratio, L/R, completely determines its values [TSH2S]. However, 
Park and Crowell 25J have observed that crystalline defects working as very small pinning sites play a crucial role 
in the vortex dynamics, namely in the critical field that resonates the gyrotropic mode. Hoffman et al 26J have 
also observed a remarkable asymmetry in the magnetization reversal mechanism which they credited to such 
small defects. Kuepper et aZ[27] have studied this mode when a (cylindrical) cavity is intentionally introduced 
in the sample. Their results clearly show how these cavities affect vortex motion as a whole, particularly its 
gyrotropic frequencies. In addition, Hoffmann et al [26] have verified that the two magnetostatic spin- wave 
modes degenerate to a single frequency when the vortex core is captured by a hole. Here, we shall show that a 
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number of facts, including some ones observed by Kuepper et al [27], may be well understood in terms of very 
basic physical properties of the vortex structure and dynamics in the presence of a cavity. Qualitatively, our 
arguments and results may be useful for other similar systems were pinning sites/defects concern. 



2 The analytical model and basic results 



Consider a magnetic dot represented by a small cylinder of radius R and thickness L (so that its aspect ratio 
L/R << 1). Assume that along the axial direction (z-axis), the magnetization M is uniform; assume also 
negligible anisotropy energy (e.g., Permalloy samples). Furthermore, if we introduce N isolated holes (each 
of them with height L and radius pn << R) in the dot, the total magnetic energy of the nanodisk can be 
approximated, in the continuum limit, bv|28j 
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where A is the exchange coupling, D is the area of the cylinder face, rh = M/M s is an unity vector describing 
magnetization along D [M s is the saturation magnetization), h m = h m (rn) = H m /M s is the demagnetizing field, 
/i ex t is an applied magnetic field (Zeeman term) and p = 1,2. The potential U, in turn, brings about the effect 
of the holes distributed throughout the nanodisk, say, ]X=i Ui(r — di) = U\(f— di) U^r — ds) ■ ■ ■ UN(r '— djy), 
with 
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Therefore, the system of a dot with N isolated antidots may be viewed as a cylinder of radius R and thickness 
L with N smaller cylindrical cavities with radius p « R, each of them centralized at di. Here, we shall study 
the case N = 1 (the treatment for N > 2 may be performed in the same way). 



Let us sketch the framework: consider a cylindrically symmetric vortex-like magnetization throughout 
the dot face, say, with the vortex core centralized at r — which is the magnetization ground state of 
a nanodisk, with L/R << 1, in the absence of both hole and external field. It is convenient to write 
rh = (sin#cos(/?, sin#sin(/3, cos#), with 6 = 6 v [r) and ip = arctan(y/x) ± n/2. The function 9 v (r) may be 
approximated by sin#„(r) = in the dot center (r = 0), while sin 8 v (r) — > 1 far away the center, | r |= r » l CK 
(Z cx = \J A/2nM s is the exchange length). In words, the magnetization consists of a small core where spins 
develop out-of-plane components for regularizing the exchange energy, and an outer region where spins are prac- 
tically confined to the dot plane face, so minimizing the (stray) magnetostatic cost. In this ground-state, the 
magnetic superficial charges in the lateral face of the dot and the magnetic volumetric charges, V ■ m, identically 
vanish yielding no contribution to h m . The cavity affects this picture as follows: it yields a less exchange energy 
for the vortex, attracting it and eventually deforming its profile, according to the potential below [T01 [TT1 [12] : 
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where s = f/R is the relative vortex center displacement, d is the distance between hole and disk centers, and 
b is a constant introduced to avoid spurious divergences whenever the vortex is centralized at the hole, i.e., 
(sR — d) = (to be estimated later) . Note that the first term inside ln-function accounts for the deformation 
of the vortex while it dislocates from the center while the remaining one takes into account the hole effect. 
For ensuring the strict validity of the above and other contributions to the total energy we must consider s 
sufficiently small, so that d and p should be also small compared to the disk radius, R. 



In addition, the distribution of magnetic charges throughout the internal edges of the hole and mainly along 
the external lateral face of the nanodisk (if the cavity is not at the center of the disk) increases the magnetostatic 
energy due to a change in the product m • n s (n s are unit vectors normal to the external lateral surface of the 
disk and internal surfaces of the cavity). Thus, the demagnetizing field, h m , can be obtained from its associated 
potential $,„ = $v + ^odgc + ^edge> m tne USU& 1 wa Y; = — V$ m . Here, $v is the magnetostatic potential 
related to the volumetric charges, while ^dgc an d ^cd^c come about from the surface charges on the external 
and internal (hole) edges, respectively. The contributions of the volumetric potential can be neglected since the 
approximations assumed above lead to V • fh = 0. Thus, the magnetostatic contribution reads [101 [TT1 [T2l [28] : 
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where F%(x) = f °° dt[J 2 (t) (e xt + xt — l)]/xt 2 . In addition, (a,f3) = (1,0) for \s ~ d/R\ < p/R (vortex center 
is inside the hole), or (0, 1) if \s — d/R\ > p/R (out from the hole). 

In the presence of an external field the following Zeeman term must be taken into account: 

V z = -7T L{R 2 - p 2 )hM 2 s + 0{s 3 ) , (5) 

with h — \H\/M S . Therefore, the effective potential experienced by the vortex reads V c s(s) — V cx + V mag + Vz- 
Depending on the relative hole- vortex centers we have the following possibilities: i) the hole is centered at 
the disk, d = (0,0), yielding a unique point that minimizes energy, s*o = (0,0); ii) the hole center is at 
d = (X, ), X > (without loss of generality) and the vortex center is far apart X, say, a = 0, so that we again 
have so = (0,0) as the unique equilibrium position (the hole practically does not affect the vortex), and; hi) 
d = X > and the vortex center is inside the hole (or very close to it), a = 1, then besides Sq it also experiences 
a new EP at s\ = (xi/R, 0), with 
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which depends on the disk and hole relative geometries and is located between the disk and hole centers. Indeed 
x\ < X by virtue of the magnetostatic effect which always tends to centralize the vortex in the disk (these 
cases are illustrated in Figure [T]). The parameter b may be estimated by equating the vortex exchange energy 
in the presence of the hole, eq. ([3]), to its normalized rest energy, Eq = 2wALhi(R/lo), in the same situation, 
-Bo hoio = 2ttAL hi(R/p), say, V ex = So hole — Eq — E COTe , where E corc — 2-kAL is the core energy, which must be 
taken into account if the vortex is out from the hole. After some algebra and assuming Iq = l ex and p > 4x, we 
get: 

InTl-^) =-4[l + ln(p/U)] , (7) 
which gives b > p (in the thermodynamical limit b = 1.04/3 |13j). 

Therefore, each non-centered defect may provide an additional EP for the vortex center, as observed in 
experiments [3 [51 [S] and predicted theoretically [101 [HI HH HE] • Namely, note that if p > p cr then the potential 
well at the hole is deeper than that at the disk center (see Fig. [TJb). In this situation, once the vortex center is 
captured by the cavity it will remain there unless a strong enough perturbation (like an external homogeneous 
field) take it to another EP. Indeed, turning h on the vortex center is shifted from its old EP, along x, by: 

Xh _ ±h(R 2 - p 2 ) 
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for h = ihy and clockwise vortex orientation. Thus, the vortex EP's with both hole and external field are 
so,h = {xh,0) and sj. j, = (xi + £/ l7 0). Furthermore, once the hole attracts the vortex to its center, there is 
an asymmetry in the vortex path as h is varied, say, from -1 to +1. At fields ~ —1 the vortex is practically 
annihilated at the disk border. As the field is gradually removed towards zero, the vortex continuously displaces 
to the disk center. However, when it gets close to the hole border it abruptly jumps inside it, at a critical field, 
/i cr i < 0, nucleating around the hole center. For a wide range of the field the vortex center remains inside the 
hole. Only at a sufficiently strong field, h CI 2 > 0, the vortex center is released from the hole, performing another 
abrupt jump (note that |/i C r2 1 > |^cri|)- This scenario is depicted in Figure [21 where the usual linear path 
is also presented for comparison. These should be compared to those obtained in the experiments of Ref.[7], 
its Figures 3 and 4 and related text; namely note that our results fit qualitatively well those reported in this work. 



In addition, there follows that the spring-like constant experienced by the vortex center reads: 
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which is clearly enhanced whenever the vortex center is inside the hole, a = 1. Results above are strictly valid 
within the rigid vortex assumption, and might be extended, for instance, to the two side charges regime. 

Besides the spring-like constant, the hole may also lead to profound modification in the vortex structure and 
dynamics. For instance, like an external field a hole may trigger the formation of non-homogeneous magnetiza- 
tion pattern at the disk and hole borders [TOl HU [T2] ; also, when its center is captured by the hole its gyrotropic 
motion is suppressed [27], Other possibilities will be discussed in what follows. It should be stressed that in the 
absence of the hole our results recover those already found in the literature [T5] . 



3 Further results and discussion 



Consider two disks, B and C, with radii Rb and Rc- Disk B has a non-centered hole, of radius pg at ds = 
(d,0), d > 0, while in disk C the hole, pc, is centered. Let A denote a disk, Ra, without hole. Their thickness 
are La, Lb, and Lc- For concreteness we also assume typical Permalloy parameters: A — 1.3x 10 -6 erg/cm = 
1.3xl0 _u J/m and M s = 800 emu/cm 3 = 800KA/m, so that the exchange length, Z cx = v / A/2nM 2 ~ 5nm. 
When vortex gyrotropic vector concerns, its gyroratio is ^/2tt — LM S /2G — 2.95 GHz/KOe = 37KHzm/A. 

Disks B and C match their counterparts 2 and 3 from Ref. [27] , with R2 = 1 /im, pi = 300 run, and d is as- 
sumed = 400 nm (not provided in Ref. [27]; we have estimated it by means of the available images); R3 = 1.5 pm 
and P3 = 300 nm, respectively. Our disk A matches their disk 1, with R± = 750 nm; all the three disks are 
L = 50 nm thick. In their work, the authors have estimated (experimental values are not presented) the gy- 
rotropic frequency of the vortex in their samples, based upon the two- vortex side charges (without hole effects), 
to be 333.5 MHz, 249 MHz, and 164.5 MHz, for disks 1, 2 and 3, respectively. Once the values above are 
believed to superestimate experimental findings in 10% — 15%, hole effects are expected to account for such 
deviations, say, in hollowed samples, like disks 2 and 3. As usual, such frequencies are given by log — k/G, 
with G = 2TrM s Lqp/"f —, where q and p account for the chirality and polarization of the vortex (the possibility 
of fractional polarization is treated below). Once the vortex is nucleated out from the hole, its gyrovector is 
completely turned on but k changes according eq. |9|), with a = 0. In this case the exchange contribution to 
k, third term in eq. ([9]), is very small compared to the magnetostatic one for these disks, so that the main hole 
contribution to the gyrofrequency comes from the factor (R 2 — p 2 )/R 2 which reads 0.91 and 0.96 for disks 2 
and 3, respectively. Thus, their respective frequencies are decreased by about 9% and 4%, in the presence of 
the hole. Such deviations appear to become larger as L/R increases as shown in Fig. [3] 

These are slight effects if compared to those in which the vortex center is captured by the hole, suppressing its 
gyrovector, once its core is removed. In this case, the vortex dynamics is deeply affected and low frequency modes 
(like the usual gyrotropic) no longer take place, as observed in experiments [27] . However, other oscillations are 
possible within the dynamical equation: 

G x x - fcx = M v x, (10) 

where x and M v = irti 2 \n(R 2 /Iq)/8AL1 2 x are the vortex position and mass (x = dx/dt, etc). This mass ex- 
pression has been adapted from Ref. [29], with Iq being the core radius. Instead of \a.(R 2 /I 2 /) a R 2 /Iq behavior 
has been recently claimed to better fitting simulations, at least for planar vortices 13 . Indeed, additionally 
to the lower gyrotropic, lug = k/G, there is an harmonic vortex oscillation around the hole with frequency 
lom — v k /M v , namely if the vortex is inside the hole. For typical samples, like those considered above, we 
find M v = 0.01 — lm e (using ln(i? 2 /Iq) and R 2 /Iq, respectively; m e is the electronic mass), so that harmonic 
oscillations take place around u>h — 10 15 — 10 16 Hz (~ 10 6 times the gyrotropic ones), which is much higher 
than currently available experimental capabilities, around some dozens of GHz }30j . 

Small vortex oscillations around a given EP naturally take place for balancing net energy of the magneti- 
zation, for instance, against an external alternating field. When this occurs some portion of the vortex central 
region can escape from the hole and possibly develops out-of-plane spin components, so developing part of the 
gyrovector. Such a gyrovector (lying, probably, at or above the GHz scale, like below) could be observed if 
a suitable alternating field were applied to resonantly turn on such a mode. In this case, additionally to the 
harmonic oscillation perpendicular to the field the vortex would also orbit around the defect (like it usually does 
around an EP). This would ensure that, whenever taken away from the hole, vortex gyrovector is developed from 
the core border towards the center, demanding gradual and continuous flipping of the normal magnetization. 
[Conversely, if no such additional orbiting motion were observed, then we should expect that vortex core (and 
gyrovector) starts to be formed from the center to the border, a fact less probable energetically, once it demands 
an abrupt 7r/2-flipping of the spin at the center]. Therefore, such experiments could also reveal the fractional 
gyrovector structure of this vortex-like configuration, provided by the dynamics of a vortex captured by a hole. 
Indeed, it should be stressed that those spin sites at the hole border (like those at the disk border) have a link 
topology different from those located at the bulk: they have only five nearest-neighbor against six from the 
bulk sites. Such a asymmetry could somewhat induce, at a tick strip around the hole border, out-of-plane spin 
component, for instance, through exchange anisotropy. In this case, a fractional gyrovector takes place, like 
below: 

M L 

G frac = — / (Vcos9 x ) x (V$)d 2 r = X G , (11) 

7 Jo 

where the integral is effectively evaluated along a small distance, S, from the hole border, r' = r — p = 0, 
towards the bulk, where VO x 7^ 0. We have taken $ = ±0 ± tt/2, <fi is the azimuthal angle on the disk plane, 



and O x is such that cos0 x — > x as r ~^ Pi \x\ < 1) an d cosO x — > for r > p + 5. A trial solution is given by 
cos(9 x ) = Tx{[{ r ^p) 2 ^S 2 ]l[{r — p) 2 + 8 2 ]} n , with n a positive real parameter (namely, n = 4 yields a smoother 
behavior for the magnetization [31]). In words, O x describes a magnetization configuration which presents a 
fractional polarization at the hole border, p x — \, falling off rapidly as we go towards disk border. Thus, 
X G (0, 1) accounts for how much the spins around the hole border are flipped towards the normal. Indeed, such 
a flipping is expected to be smaller as the hole increases, so that above a critical value no net magnetization 
perpendicular to the disk plane takes place anywhere. In the case of a static vortex this critical hole size goes 
around ~ 0.3Z ex [32] . but could be largely increased in the dynamical case, once the vortex center can now moves 
towards hole border and even escape outside. Therefore, whether this fractional polarization actual takes place 
seems to strongly depend on the relative location of the vortex center to the hole border at a given time, varying 
according the distance between them changes. However, if the hole is very large the displacement of the vortex 
center to the hole border may be highly energy costing leading eventually to a deep deformation in its profile 
or even its annhilation, so that its central region is kept inside the hole yielding no net polarization at all. This 
seems to be the case of sample 3 studied in Ref. [27] , whose observations have led the authors to conclude that 
the gyromode is completed suppressed, once the vortex core is captured by the centralized hole. Nevertheless, 
just for comparison, consider that the gyrovector was not completely turned off, but was reduced to, say, 10% 
(X — 0.1, i.e., maximum spin flip, O x | max = 9°). In this case, sample with p = 300 nm and R = 1,5pm, the 
frequency associated to the fractional gyrovector reads 

OW = fe/G frac ~lQGHz, (12) 

increasing as x is diminished. Even in this situation of fractional gyrovector our analysis may explain the reason 
why in the work of Ref. [27] no trace of the subGHz gyrotropic mode was observed: the applied field frequency 
is far below that required for resonantly excite the fractional mode, as eq. ([12]) predicts. 



4 Conclusions and Prospects 



We have realized that the introduction of a hole into a ferromagnetic nanodisk, with small aspect ratio, may 
deeply modify the structure and dynamics of a vortex-like magnetization. Such effects come about once a hole 
change the effective potential experienced by the vortex in such a way that the vortex center is attracted towards 
the hole. As a vortex is captured by a hole its out-of-plane components, accounting for its polarization, are 
greatly diminished or even vanished. Such a capture may be clearly viewed in the plot of the vortex center 
position against a uniformly varying external field, Figure [21 A moving vortex whose center is not inside a hole 
experiences a decreasing in its associated gyrotropic frequency, which may become pronounced as the disk and 
hole sizes are comparable. 

We have also raised the interesting possibility of the fractional polarization of moving vortex. This would 
happens if the vortex center were inside the hole, but sufficiently close to its border, so that some portion of 
the out-of-plane magnetization show up. Among other issues, its observation could answer how the vortex core 
is formed as it is released from a hole. Our results could have some relevance to the study of structural defects 
in nanomagnets, once a analytical model for understanding such structures in these frameworks is still lacking. 

A very interesting result to be addressed is the analysis of the vortex motion in the presence of the hole and/or 
external field. Analytically, this task demands the resolution of complicated differential equations obtained from 
eq.([ID]) with the force, F = M„x, given by F = — VT4ff- Simulations have revealed the richness of such dynam- 
ics, including the possibility of vortex core reversal (switching) triggered by vortex-hole interaction 33 . 
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5 Figure Captions 



Figure 1: [Color online] Typical plots of the potential V c s/A as a function of the vortex position, s. Here, 
we have taken R = 1000 nm, L = 50 nm (hole sizes are shown in the figures), a) centralized defects produce 
deeper and sharper potential wells as their sizes increase. Here, the graph corresponds to the normalized 
potential so that the bottons are depicted at the same normalized value, b) each eccentric hole induces a 
new minimum for the potential, which may be the absolute if the hole size is large enough. Since the disk 
center always attracts the vortex, the minimum due to an eccentric defect does take place at a given point 
between the disk and hole centers (here, the holes are centralized at x = X/R = 0.4). 

Figure 2: [Color online] The vortex equilibrium position as function of an external varying field. The linear 
behavior is observed without the hole. A centralized cavity deeply changes vortex dynamics as the field 
is continuously increased from -1 to +1, namely there appear abrupt jumps when its center is captured 
(released) by (from) the hole (blue dotted lines). Note also that the graph slop is greatly lowered whenever 
the vortex center is inside the defect (red line), exhibiting the stronger effects of the net potential (for 
improving visualization, we have shown only the central region of the graph, where the differences take 
place) . Compare with related results from Ref. [7] . 

Figure 3: [Color online] How the gyrotropic frequency of the vortex, u>q, is influenced by the presence of 
a centralized hole (p = 100 nm). The change is slight for very large disks, but becomes pronounced as the 
disk radius is decreased, namely when disk and hole sizes becomes comparable. 
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